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昆虫 抗 药性 的 治理 策略 : 
数学 模型 的 提出 
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WE ”本文 用 计算 机 模型 分 析 及 评价 了 昆虫 抗 药 性 的 三 种 治理 策略 ; 顺序 纶 用 ，“ 高 杀 死 ”策略 及 棋盘 
式 用 药 。 假 设 昆 虫 具有 一 定 的 生物 学 特性 ; 如 抗 福 为 单 基因 >, 单 抗 性 , 抗 性 基因 为 半 显 性 ,两 福生 殖 , 虹 虫 有 
扩散 习性 * 感 往 系 在 不 施 药 情况 下 有 微细 的 适应 优势 等 。 分 析 了 用 不 同 痢 量 造成 的 不 同 的 死亡 率 : 不 周 的 反 
选择 作用 * 不 同 的 基因 稀释 作用 等 情况 对 三 个 策略 的 影响 。 结 果 明 示 出 > 在 这 些 人 情况 下 ,用 3 一 5 种 杀 虫 剂 隔 
代 涛 序 轮 用 加 上 微弱 的 反选 择 作 用 是 最 为 有 效 的。 “高 杀 死 策略 加 上 一 定 的 稀释 作用 也 极为 有 效 。 棋盘 式 
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昆虫 抗 药性 治理 在 近 几 年 来 已 逐渐 成 为 引 人 注 目的 问题 。 在 1975 FRA Te Bt 
究 解毒 酶 ,企图 从 它 的 解毒 机 制 出 发 来 克服 昆虫 抗 药 性 。 事 实证 明 ， 这 种 方法 并 不 成 功 ， 
搞 药 性 一 旦 形成 后 ,往往 就 难以 消除 。 因此 上 且 前 对 抗 性 的 研究 正 转 向 阻止 或 延缓 抗 性 的 
形成 。 这 就 是 抗 性 治理 的 研究 目标。 - 
昆虫 抗 药性 的 形成 ,实质 上 是 昆虫 群体 中 抗 性 基因 频率 增加 的 过 程 。 增加 的 原因 主 
要 是 杀 虫 齐 的 选择 作用 。 抗 任 治理 的 目的 ,就 是 设法 不 使 抗 性 基因 增加 ,或 延缓 其 增加 的 
速度 ,同时 又 要 控制 害虫 为 害 ( 张 宗 炳 ,，1986、1987a; Bishop, 1982; Conway & Comins, 
1979; Keiding 1977; Muggleton, 1982). 当然 ， 不 用 化 党 药 州 是 最 彻底 的 抗 性 治理 手 
艘 , 但 这 在 目前 还 难于 做 到 。 抗 性 治理 的 基本 原理 有 三 种 : (ORERE; CORALE 
择 作 用 , 即 不 施 药 时 感性 昆虫 更 强 的 适应 性 ;(3) 利 用 稀释 作用 , 即 让 防治 区 周围 的 感性 害 
虫 迁 人 并 与 抗 柱 个 体 杂 交 。 在 不 同 的 条 件 下 ,这 三 者 的 作用 强 弱 是 不 同 的 ,应 综合 郑 虑 害 
虫 生物 学 特性 和 具体 条 件 ,才能 收 到 最 佳 的 效果 ( 张 宗师 ,1987a、1987b) ,本 文 就 是 利用 种 
. 群 遗 传 学 原理 ,采用 计算 机 模拟 方式 ,定量 地 对 几 种 抗 性 治理 策略 进行 了 评价 。 由 于 害 出 
-生物 学 特性 对 抗 性 形成 有 明显 影响 ( 张 宗 炳 ，1987b; Tabashnik & Croft, 1982), 模拟 时 
对 这 一 点 必须 加 以 其 虑 ,本 文中 目标 昆虫 共有 以 下 的 生物 学 特性 : 〈1) 抗 性 是 由 单 基因 控 
; 制 的 ,并 且 是 单 抗 性 。(2) 该 基因 是 半 显 性 遗传 ， 抗 福 水 平 高 。(3) 害 虫 是 两 性 生殖 。(4) 
- 害虫 有 较 强 扩散 能 力 。(5) 害 虫 无 世代 重 迭 现象 。(6) 感 性 害虫 在 不 施 药 时 比 抗 性 害虫 有 
微弱 的 优势 , 即 存 在 反选 择 压 。(7) 初 始 抗 性 基因 频率 低 。 
| 我 们 使 用 的 3 一 5 种 杀 虫 剂 ,没有 交互 抗 性 ,但 也 不 要 求 有 负 交 互 抗 性 (例如 ， 可 选用 
一 种 有 机 磷 , 一 种 拟 除虫菊 栈 ， RANE » DKS — EE PRIS Jo 它们 的 抗 性 基因 不 能 


本文 于 1987 年 6 月 收 到 。 
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”有 连锁 遗传 的 现象 。 

根据 以 上 的 基本 假设 ,我 们 模型 中 主要 包括 两 组 公式 : ”第 一 组 公式 描述 种 群 数量 的 
变化 , 它 的 核心 是 以 togistic 晶 线 型 的 差分 方程 来 描述 种 群 的 增长 ,再 进一步 计算 施 药 后 
不 同 基因 型 个 体 的 存活 率 。 第 二 组 公式 则 是 在 此 基础 上 计算 出 抗 性 基因 频率 的 变化 。 主 
要 公式 为 : 

第 * 代 种 群 的 初生 数量 : 

NG —1,2)+ BR 
NG, 1) = OCEN 《1 
K 


其 中 NG 一 1,2) He 一 1 代 成 熟 个 体 数 ; BR 为 繁殖 力 ， 即 没有 密度 制约 情况 下 相 邻 世 
代 成 虫 数 量 之 比 。 它 实际 上 是 昆虫 产 卵 县 与 自然 成 活 率 的 乘积 ， 本 模型 中 取 3.5; 玉 为 环 
境 容量 ,本 模型 中 取 为 10000。 本 公式 实际 上 是 常用 的 logistic 曲线 的 差分 形式 。 

施 药 后 昆虫 的 存活 率 由 施 药 剂量 决定 。 为 了 使 死亡 率 与 药物 剂量 的 对 数 成 为 线性 关 
系 , 我 们 对 死亡 率 进行 了 机 率 值 变 换 。 详 细 计 算 过 程 请 参见 昆虫 毒 理 学 。 

由 于 存在 反选 择 压 ,我 们 认为 杂 合 子 的 存活 滨 比 感性 纯 合 子 要 低 一 些 , 应 再 科 以 系数 
1-AS, KEM AS 为 反选 择 压 。 至 于 施 药 后 的 存活 率 , 假 设 不 同 基 因 型 其 死亡 机 率 值 与 
剂量 对 数 的 回归 线 的 斜率 相同 ,只 是 截 距 不 同 ,因此 我 们 用 P-H 来 表示 杂 合 子 (SR) 的 
死亡 机 率 值 ,其 中 PP 为 相同 剂量 下 感性 纯 合 子 《55) 的 死亡 机 率 值 ,日 为 杂 合子 (SR) it 
性 显著 度 , 本 文 取 为 2.33, 这 大 约 相当 于 感性 纯 合 子 《55) 接近 全 部 死亡 时 , 杂 合 子 (SR) 
死亡 一 半 。 

对 于 抗 性 纯 合 子 (RR)， 它 的 存活 率 也 应 乘 以 系数 1-AS， 这 意味 着 抗 性 纯 合 子 
(RR) 与 杂 合 子 (SR) 生活 力 相同 。 它 的 施 药 后 存活 率 为 1, 即 在 实际 使 用 的 药物 浓度 下 
它 不 被 杀 死 。 

迁 人 量 可 有 几 种 不 同 的 算法 ; 但 如 果 我 们 假设 周围 环境 中 该 种 昆虫 种 群 数量 相对 稳 
定 的 话 , 则 认为 迁 人 量 为 固定 常数 是 合理 的 。 

假设 迁 人 昆虫 的 抗 性 基因 频率 为 0.01 固定 不 变 , 则 不 难得 出 上 世代 种 群 中 抗 性 基因 
频率 为 : 
2M(G,1) 十 M(152) + 0.021 (2) 

2N(t,2) 
其 中 M (451) 为 :世代 成 熟 种 群 中 抗 性 基因 纯 合子 个 体 数 ; M(1,2) 为 杂 合子 个 体 数 : 
为 迁 入 量 ; NCr ,2) 为 + 世代 成 熟 期 种 群 总 个 体 数 (包括 迁 人 )。 

模拟 过 程 中 ， 模 型 输入 的 固定 参数 有 :繁殖 力 BR 一 3.5; 迁 人 种 群 抗 性 基因 频率 
0.01。 可 变 参 数 有 ， 反 选择 压 AS, EAE !， 药 物 对 种 群 的 杀 死 率 CD。 由 于 种 群 抗 
性 水 平 不 断 变化 ,为 了 维持 相对 稳定 的 CD 模型 中 要 求 它 的 波动 不 超过 5 %) ,用 药剂 量 
也 是 变化 的 。 但 抗 性 水 平 很 高 时 ,指定 的 CD 可 能 达 不 到 (因为 抗 性 纯 合子 不 会 被 杀 死 )， 
显然 这 样 的 假设 更 接近 实际 使 用 情况 。 模型 的 初 值 一 般 为 种 群 初始 数量 N(0,2) 一 
1000， 这 也 是 我 们 选 定 的 经 济 阅 值 ; 抗 性 基因 初始 频率 PCO) = 0.01。 

模型 输出 采用 图 形 的 方式 。 图 中 纵 坐 标 为 抗 性 基因 频率 P 和 种 群 数量 No RHA 
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世代 数 。 每 组 的 2 一 3 个 图 坐标 相同 。 . . 

图 1 为 抗 性 增长 的 模拟 。 其 中 (Ca) 为 迁 人 的 影响 。 显 然 迁 人 的 增加 只 能 稍稍 推迟 
抗 性 增长 。 虽然 药物 浓度 不 断 提 高 ， 但 能 控制 种 群 数量 不 超过 经 济 阔 值 的 时 间 只 有 5 一 
10 Ro (b) 为 反选 择 压 的 影响 。 反 选择 压 的 增加 对 抗 性 增长 速度 几乎 没有 什么 作用 ,只 
是 使 种 群 最 高 数量 有 所 降低 ,但 这 没有 意义 ;因为 它 早 已 超过 了 经 济 阐 值 。Cc) 是 选择 压 
的 影响 ,这 里 可 以 看 到 一 个 有 趣 现 象 , 即 从 开始 就 提高 药剂 浓度 后 ， 可 使 抗 性 维持 在 很 低 
水 平 上 。 这 觉 明 在 有 迁 人 的 情况 下 高 杀 死 策略 很 可 能 是 有 效 的 。 关 于 这 一 点 我 们 后 面 还 
村 详细 讨论 《Wood & Mani，1981: Taylor & Georghion, 1979)o 
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图 1 抗 性 增长 每 代 施 药 
(a) 迁 人 的 影响 :43 一 0.05 CD= 0.8, (1) 7 一 6 (2) L= 500 
(b) 反选 择 压 的 影响 ， CD = 0.8 1= 100, (1) AS=0.05 (2) AS=0.2 
(c) EBM: [= 100 AS = 0.05, (1) CD= 0.8 (2) CD = 0.95 


从 抗 性 增长 模拟 可 见 , 一 般 情 况 下 抗 性 增长 是 很 快 的 ,第 2 一 3 代 就 会 超过 0.2， 能 够 
控制 种 群 数 量 的 时 间 不 超过 10 代 。 因此 仔细 研究 抗 镍 治理 策略 防止 这 种 现象 发 生 是 十 
分 必要 的 。 
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(a) (b) 
图 2 抗 性 不 退 模拟 ,不 施 药 , 初 值 PCO) = 0.2 


(a) FAME: AS = 0.05, (1)1=0 (2)1=100 (3) 1= 500 
Cb) 反选 择 压 的 影响 : I= 100, (1) AS = 0.05 (2) AS=0.1 (3) AS= 0.2 
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AZAR BURN, TARR RES To 抗 性 基因 频率 从 0.2 pag 
0.01 Æ AS= 0.05, (== 0 的 情况 下 要 60 多 代 , — SOOM HE 40 SR, 在 = 100, 
AS = 0.2 的 情况 下 也 要 15 RAH 

3A RAR RNR, PARR 是 指 采 用 高 药剂 量 ， 以 求 把 杂 合 子 个 体 
全 部 杀 死 ,这 样 就 减少 了 抗 性 基因 的 数量 。 在 有 一 定量 迁 和 人 ,感性 基因 稀释 的 情况 下 ， 就 
可 能 阻止 种 群 抗 性 的 上 升 。(a) 是 迁 人 量 的 影响 。 由 于 杂 合 子 全 被 杀 死 ， 所 剩 下 的 抗 性 
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(a) i (b) 

图 3 高 杀 死 策略 ，A5 = 0.05, SSH 为 如 合子 存活 率 

(a) TAME, SrH = 0, (1) om 3 (2) 710 

(b>) 杂 合 子 存 活 率 SrH 的 影响 : I 100, (1) SrH = 0.22 (2) SrH = 0.21 
纯 合子 极 少 ( 在 所 给 的 初 值 下 只 有 0.1 头 ), 迁 人 量 影响 不 大 。(b) ERAT HER (StH) 
的 影响 。 模 型 对 这 个 参数 极为 敏感 , 当 存活 率 从 0.22 变 到 0.21 时 ， 模 型 的 行为 有 了 根本 
的 改变 ， 从 不 能 控制 抗 性 和 数量 变 为 能 够 控制 。 模 型 运行 说 明 , 只 要 有 一 定 迁 入 Mitt 
治理 的 观点 看 ， 高 杀 死 "策略 是 有 效 的 ;但 从 环境 保护 的 观点 看 ， 高 剂量 使 用 杀 虫 剂 记 带 
来 的 副作用 也 许 是 不 能 接受 的 。 因 此 它 的 使 用 就 受到 了 很 大 限制 。 
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图 4 棋盘 式 用 药 ( 两 种 药 ) 
(a) 迁 人 量 的 影响 ， 4 一 0.05 CD=0.8, (1) 7 二 100 (2) P= 5 的 
(b) 杀 死 率 的 影响 : A43 = 二 0.05 7 =100, (1) CD=0.8 (2) CD 一 0.9 
(c) 反选 择 压 的 影响 : CD 二 0.8 £5100, (1) AS=0.05 (2) AS = 9,1 
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图 4 是 对 “棋盘 式 ” 用 药 的 模拟 。 这 种 策略 是 象 国际 象棋 的 棋盘 一 样 , 黑 格 用 一 种 药 ， 
” 色 格 用 另 一 种 药 。 这 桩 不 同 的 抗 性 基因 可 以 互相 稀释 。(a) 图 是 迁 人 量 的 影响 。 可 以 看 
到 [一 100 时 仍 不 能 阻止 抗 性 形成 ,而 当 1 一 500 时 可 以 阻止 。 但 这 相当 经 济 阔 值 的 
一 半 , SLA EEDA. 而 且 平衡 时 的 抗 狂 基因 频率 为 0.2 左右 ， 也 太 高 了 一 
3, O) 是 施 药 剂量 的 影响 。 如 果 药 物 杀 死 率 CD 一 0.9， 则 可 以 把 抗 性 和 数量 都 控制 
在 低 水 平 上 。 这 里 又 显 出 了 “高 杀 死 ”策略 的 效果 。(c) 是 反选 择 压 的 影响 。 当 4$ 一 0.1 
时 ,种群 数 量 NN 可 以 控制 ,而 抗 性 基因 闫 率 P 则 波动 于 0.3 一 0.4 之 闻 ， 这 种 波动 是 由 于 施 
药 浓 度 变 化 引起 的 。 在 本 图 中 ,三 个 对 照 (1) 的 条 件 是 相同 的 。 
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”图 5 五 种 药物 轮 用 ,每 种 用 一 代 
(a) ZARE: 443 一 0.05 CD=0.8, (1) 7 一 0 (2) 7 一 2 
(b) 杀 死 素 的 影响 : AS= 0.05 1 二 100, (1) CD 一 0.5 (2) CD=0.7 
(O 反选 择 压 的 影响 ， CD 一 0.51= 100, (1) AS= 0.05 (2) 4S =0.1 


图 5 是 对 “顺序 轮 用 ”的 模拟 。 在 本 模拟 中 采用 5 种 没有 交互 抗 性 的 药物 ,每 种 用 一 
代 。(a) 是 迁 人 量 的 影响 。7 一 0 时 虽然 最 终 仍 不 能 阻止 抗 竹 上 升 ， 但 已 大 大 延长 了 控 
人 币 数 量 的 时 间 ( 可 达 25 代 )。 而 只 要 有 一 点 稀释 作用 ,即使 仅 为 1 一 2， 也 可 使 抗 性 和 数 
量 都 维持 低 水 平 (此 时 图 中 的 N 与 x 轴 几 乎 重合 )。(Cb) 为 杀 死 率 的 影响 。 由 于 此 时 还 
ABS (1 一 100)， 仅 当 CD 一 0.5 时 才 会 引起 抗 性 的 上 升 。CD 一 0.7 时 已 能 控制 搞 
性 ,但 基因 频率 波动 仍 较 大 。(c》 为 反选 返 压 的 影响 。AS 一 0.1 时 抗 性 虽 仍 在 上 升 ， 但 
已 极 慢 ,100 代 后 仅 为 0.1 左右 。 种 群 数量 较 多 ,是 杀 死 率 太 低 的 原因 (CD 一 0.5)。 

图 5 是 轮 用 中 使 用 代数 的 影响 。 显 然 连用 代数 越 少 越 好 。 在 (0) 中 可 以 看 到 一 个 
有 趣 的 现象 : 此 时 抗 性 基因 频率 上 升 和 下 降 的 速度 凑 不 多 。 这 启示 我 们 ， 当 害虫 已 对 某 
种 农药 发 生 较 强 抗 性 时 ,只 要 换 用 另 一 种 药 把 种 群 数量 压低 ,再 引 人 相 当 数 量 的 感性 个 休 
进行 稀释 , 则 可 在 5 一 7 代 内 使 抗 上 基因 频率 从 0.6 一 0.8 降 到 0.01 左右 。 
图? 是 三 种 策略 的 比较 。 为 了 增加 可 比 性 ,本 次 模拟 中 改 为 用 药剂 量 固定 ,这样 各 基 
因 型 存活 率 固定 ,但 整个 种 群 的 杀 死 率 “CD 不 再 是 固定 的 了 。 (C) 为 高 杀 死 策略 , 它 的 
效果 很 好 ,但 从 防止 污染 的 角度 看 很 可 能 是 不 可 取 的 。(Cb) 是 棋盘 式 用 药 , 它 不 能 达到 抗 
”性 治理 的 目的 o 而 在 同样 的 条 件 下 ， 旺 序 轮 用 (c》 却 可 以 得 到 较 好 效果 : ”数量 稳定 在 
200 KEF, RAMEE 0.1 以 下 波动 。 因此 对 具有 本 文 所 假设 的 生物 学 特性 的 昆虫 来 
说 ,选择 3 一 5 种 没有 交互 抗 性 的 农药 进行 顺序 轮 用 是 一 种 较为 理想 的 抗 性 治理 策略 。 
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图 7 三 种 策略 的 比较 (= 100 43 一 0.05 药物 剂量 固定 


C) 高 杀 死 策略 仅 用 一 种 药 SrS (感性 个 体 存活 率 ) 一 0.01，S5H《〈 杂 合 了 于 个体 存活 率 ) 0.1 
(b) 三 种 药 , 棋 盘 式 使 用 ，SrS = 0.1 SrH = 0.85 
(co) 三 种 药 , 顺 序 轮 用 ，SrS 一 0,1 StH = 0.85 


本 文昌 仅 从 理论 上 进行 了 计算 ,但 它 所 依据 的 是 已 经 过 验证 ,获得 了 公认 的 群体 遗传 
学 原理 ,因此 是 较 可 靠 的 《Bishop，1982)， 可 以 对 害虫 防治 工作 起 一 定 指导 作用 。 对 于 
具有 其 他 类 型 生物 学 特性 的 害虫 我 们 也 将 进行 类 似 工 作 ， 并 已 开始 对 其 中 的 一 些 重要 结 
论 进行 实验 验证 。 结 果 将 另 文 报道 。 
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ANALYSIS AND SIMULATION OF THE STRATEGIES OF !NSECTICIDE 
RESISTANCE MANAGEMENT WITH A MATHEMATICAL MODEL 


S. G. Li J3 T. CHANG J. M. Yan 


(Depariment of Biology, Peking University, Beijing) 


A mathematical model is set up to analyze and simulate three different strategies of in- 
secticide resistance management: the sequential use of 3—5 insecticides; the “high kill” policy, 
with immigration of susceptibles; and the “patch” or “checker-board” system of insecticide ap- 
plication, The results show that the sequential use of 3—5 insecticides, each for one genera- 
tion, gives the best result provided that there is no cross resistance among the insecticides. 
“High kill” policy with immigration of certain mumber of susceptibles is also an effective 
means, though it may not be acceptable from the viewpoint of integrated pest management. 
“Patch” system of insecticide application is generally a failure, unless the immigration rate 
is very high (about one half of the initial population). 
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